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We investigate the Josephson effect in TNT and NTN junctions, consisting of topological (T) 
and normal (N) phases of semiconductor-superconductor ID heterostructures in the presence of a 
Zeeman field. A key feature of our setup is that, in addition to the variation of the phase of the 
superconducting order parameter, we allow the orientation of the magnetic field to change along 
the junction. We find a novel magnetic contribution to the Majorana Josephson coupling that 
permits the Josephson current to be tuned by changing the orientation of the magnetic field along 
the junction. We also predict that a spin current can be generated by a finite superconducting phase 
difference, rendering these materials potential candidates for spintronic applications. Finally, this 
new type of coupling not only constitutes a unique fingerprint for the existence of Majorana bound 
states but also provides an alternative pathway for manipulating and braiding topological qubits in 
networks of wires. 
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I. INTRODUCTION 

The recent indications in experiments with 
semiconductor-superconductor ID heterostructures 
"J 
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of the existence of Majorana bound states (MBS) 
23| have also intensified the pursuit of designing a 
quantum c omp uter based on topologically protected 
qubits [13, l25| . The existence of a pair of spatially 
separated MBS in these systems will signal the ap- 
pearance of a topological qubit. Here, the two- level 
system arises from the emergence of a zero-energy 
quasiparticle leading to a doubly degenerate ground 
state. The topological character of the qubit reflects an 
intrinsic particlc-holc symmetry of the specific system, 
which provides the protection against external sources 
of decoherence and noise as long as neither a gap closes 
nor quasiparticles are excited [26l43l| . 

The potential applications of this class of topologi- 
cal qubits require, on one hand, finding unique finger- 
prints that unambiguously confirm their existence and, 
on the other hand, developing techniques that permit 
their manipulation with an eye to quantum informa- 
tion processing (32|-|43j. The majority of previous pro- 
posals for their detection p4Tj52l | rely on tunneling and 
transport features of the zero-energy fermionic excita- 
tion. More recently, alternative experimental routes have 
been explored, based on unconventional Josephson sig- 
natures [53l | . distinct from the well-known 47r-periodic 
Josephson effect typical for MBS 

In this manuscript, we report on new results concern- 
ing the Josephson effect in double-junctions of the type 
TNT and NTN, consisting of topological (T) and normal 



Zeeman field 




FIG. 1: (Color online) Platform for realizing topological 
qubits. A ID semiconducting quantum wire with strong spin- 
orbit coupling is placed on top of a bulk s-wave superconduc- 
tor. Proximity effects induce superconducting pairing to the 
quantum wire. Majorana bound states appear when we addi- 
tionally apply a Zeeman field. The junctions that we consider 
should be fabricated by this type of hybrid structures. 



(N) phases. Each of the segments of the junction are 
constructed from heterostructures of ID quantum wires 
(e.g. InSb) deposited on top of a bulk s-wave super- 
conductor (e.g. Nb), in the simultaneous presence of a 
Zeeman magnetic field (Figure [lj. We demonstrate that 
there exists a contribution to the Josephson coupling of 
these systems which can be tuned by varying the orienta- 
tion of the Zeeman magnetic field along the junction. For 
two neighbouring MBS, this novel "magnetic" Josephson 
coupling in its simplest version is proportional to 



H M " OC ljlj r COS 



if I - lf r 



(1) 
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where 7;. r correspond to the Majorana operators satisfy- 
ing {7i, r ,7!,r} = 1 and {7z,7 r } = 0, correspond to 
the angles defining the orientation of the Zeeman field 
and (pi,r represent the superconducting order parameter 
phases that the MBS see. 
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We observe a 47r-periodicity also for the Zeeman field 
phase difference i?/ — i? r , that could provide a character- 
istic signature of the interacting MBS. The recent ex- 
perimental results 0] supporting the observation of the 
MBS Joscphson effect in these heterostructures, suggest 
that our setup constitutes a realistic, unambiguous and 
accessible way to detect the possible existence of MBS. 
Moreover, this extra term may offer alternative routes for 
braiding and manipulating MBS in networks of quantum 
wires, as for example in Y-junctions [38l l4lj. by adiaba- 
tically varying the relative phases of the Zeeman field. 

II. MODEL HAMILTONIAN AND BULK 
SINGLE PARTICLE EIGENFUNCTIONS 

In order to study the Joscphson effect of the junctions, 
we model the semiconductor-superconductor ID hetero- 
structure with the following Hamiltonian 

71 = \ / dX ^( X W(P, x )$( x ) > ( 2 ) 

where we have introduced the Bogoliubov - de Gcnncs 
Hamiltonian density 

H(p,x) = vp<j z + \B(x)\e- M ^ T ^T z a x 

- \A(x)\e-^^T y a v +^-^r z , (3) 

and the enlarged Nambu spinor 

*t(a:) = ( rl>l(x) ^(x) ^(x) ) . (4) 

The factor 1/2 in front of the Hamiltonian in Eq.([2]) ac- 
counts for the doubling of the degrees of freedom when 
we introduce the extended Nambu spinor, which is com- 
posed by electronic creation and annihilation operators 
ip\(x) and ip a {x) of spin projection a =^, I- To represent 
the possible terms of the Hamiltonian within this Nambu 
formalism, we have made use of Kronecker products of 
the Nambu-spacc r and spin-space er Pauli matrices. The 
velocity v corresponds to the strength of the spin-orbit 
interaction that is oriented along the z— axis, B{x) repre- 
sents the Zeeman magnetic field that is allowed to rotate 
freely in the x — y plane, A(.t) = (Asr(x), Aq (x)) defines 
the proximity induced superconducting order parameter 
of the quantum wire written in a convenient vectorial 
form and finally the last term corresponds to the kinetic 
energy of the electrons measured from the chemical po- 
tential fi. 

For a bulk system, B(x) and A(x) become homoge- 
neous and we can Fourier transform Eq.([3]), obtaining 

U{k) = vhka z + \B\e- li>T ^^T z a x -\A\e~ lVT 'T y (Ty 

+ Ur-^K (5) 



Correspondingly, the extended Nambu spinor becomes 

¥ k = ( Vlt.^L ^-fct ^-M )■ Thc B °g° liubov oper- 
ators that diagonalize the above Hamiltonian, are gen- 
erally of the form 7fc = u fef V> fet + it^ + V^->=T + 
v k \i}y_ k \- Majorana operators, satisfying j k = jj., may 
only occur for inversion symmetric momentum space 
points (fc = —k) for which we may have linear combina- 
tions of the type ip k<7 ± ^_ ka = ^ ka ± with a =t, |. 
In the case of a bulk system this can take place only for 
k = 0. In fact, the parameter regime where the energy 
spectrum E(k) shows zeroes for k = 0, indicates a phase 
boundary between the topological trivial and non-trivial 
phases. For thc Hamiltonian of Eq.©, one directly ob- 
tains that for \B\ > \/\A\ 2 + the system is in the 
topological phase with one MBS for k = 0, while for 
| £5 1 < y/\ A| 2 + /i 2 the system is in the topologically triv- 
ial phase with a zero or an even (in general) number of 

mbs dum. 

In our case we set the chemical potential equal to zero 
/j = 0, which permits us to consider a truncated version 
of the Hamiltonian of Eq.([5]). Specifically, we drop the 
kinetic energy term and our truncated Hamiltonian reads 



H tr (k) = vhka z + \B\e-^ T ^T z c7 x - \A\ e -^T y a y . (6) 



Considering the latter Hamiltonian in order to study the 
Josephson effect is naturally justified, since we want to 
examine topological properties that are related to the 
k = point. The truncated Hamiltonian is the linearized 
version of our initial one, about k — 0. Nevertheless, 
as already discussed in Ref. [53| the models of (|5|) and 
([5]) show a discrepancy concerning the characterization 
of the topological and normal phases. Specifically, in 
the linearized model the topological phase occurs when 
|A| > \B\ and the normal when |A| < \B\. 

For studying the junctions we shall consider that the 
superconducting order parameter and the magnetic field 
remain constant for each segment, having the follow- 
ing spatial profile B(x) = Bi + (B m — Bi) <d(x — x a ) + 
{B r - B m ) &{x-x b ) and A{x) = A z + (A m - A z ) Q(x- 
x a ) + (A r — A m ) <d(x — Xb), where 6(x) is the Heaviside 
function and the label s = l,m,r denotes the left, mid- 
dle and right segments. Every segment can be described 
by the Hamiltonian of Eq.©. Therefore it is eligible to 
determine the bound state eigenfunctions for the latter 
Hamiltonian. Since we are looking for bound states, we 
set k = ik. Thc 4-component single particle bulk wave- 
functions ^(x) = e KXl $>(K,) are readily obtained and read 
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\ Kl (E);B,A) 
\k 2 {E);B,A) 
\k 3 (E);B,A) 
\k 4 (E);B,A) 



1 

^= (*cos(^) e" 



= ±( 



sin 



e '~ 



-« cos 



iJ£±l 



,- <p— a 



(7) 



(^)e+^, +sin(^)e+^) T , (8) 



s_, sm {—) e 



-l COS I 



je+'V,] , (9) 



-I sm ( — ) 2 g_ 



, (io) 



where we have made use of the definitions uj± = ia± + 
tt/2, tanha± = E/\\A\ ± and s_ = sign(|A| - 
The above wave-functions are characterized by the cor- 
responding "wave- vectors" Ki(E) = +n + (E), K2(E) = 
-k+(E), k 3 (E) = +K-(E) and k a (E) = -K-(E) where 

K± (E) = yJ(\A\±\B\) 2 -E*/vh . (11) 

For every bulk segment s = l,m,r the total wavefunction 
corresponding to energy E, retains the form 

4 

tff (x) = c s , n {E)e K '^ E > \K.. n {E); B S ,A S ) , (12) 

n=l 

where the coefficients c s;n (i?) need to be determined by 
imposing appropriate matching conditions at the two in- 
terfaces where the three segments meet pairwise. By ta- 
king into account that the junction extends to infinity 
on the left and right parts, we conclude that the total 
wave-function of the system can be written as 

*f(ar) = Cl ,i(E)e +K ^ x -^ \+k+(E); B u A,) 

+ c 3 ,i(E)e +K -^ x -^ \+k-{E);B,,Ai), (13) 

**(x) = c 2 , m (E)e- K +^ x - x ^\-K + (E);B m ,A m ) 
+ c 4 , m (E)e- K -^ x - x ^ \- K .(E);B m ,A m ) 
+ Cl . m (E)e +K ^ E ^-^ \+K + (E);B m ,A m ) 
+ c 3 , m (E)e +K -^ x - x ^ \+k-(E); B m , A rn ) ,(14) 

= c 2 AE)e-^ E K x - x ^ \-K+(E);B r ,A r ) 

+ c 4>r (E)e- K -^ x - x ^ |-K_(£);B r ,A r ). (15) 

The rest of the appearing coefficients in the above ex- 
pressions will be determined by demanding continuity for 
the wave-function along the junction. There is one con- 
tinuity equation for each interface point and they read 
tff (x a ) = V E (x a ) and V E (x b ) = V E (x b ). The aris- 
ing equations define a homogeneous system of equations 
from which we may retrieve the coefficients c Sin (E) and 
additionaly the energy eigenvalues of the system. 



III. TOPOLOGICAL - NORMAL - 
TOPOLOGICAL (TNT) JUNCTION 

We now proceed with examining specific junction se- 
tups. In this paragraph we shall focus on interfaces 
with a topological-normal-topological sequence of phases. 
We may recall from the previous section that a segment 
s = l,m,r is in the topological phase if |A S | > \B S \, 
while the opposite condition holds for the normal phase. 
In order to obtain a double junction with the desired 
phase sequence per segment, we shall consider the follow- 
ing profile for the magnetic field and the superconduct- 
ing order parameter \Bi\ = \B\, \B m \ = |A|, \B r \ = \B\, 
\Ai\ = |A|, A m = \B\ and |A r | = |A|. As a matter 
of fact, |A| and \B\ constitute the values of the order 
parameter and the magnetic field in the topological seg- 
ments. To make a connection to the results that we have 
obtained so far, the specific choice of values, renders only 
the variable s_ = sign(|A| — \B\) spatially dependent, 
which enters in Eq.© and Eq. (fi"0)) . 

After applying the matching conditions for the conti- 
nuity of the wave-functions along the junction, we may 
obtain the bound states wave-functions and the energy 
eigenvalues of the system. Due to the particle-hole sym- 
metry of the Bogoliubov - de Gennes Hamiltonian, we 
expect for every positive energy eigenvalue, an accompa- 
nying negative one. For a finite distance d between the 
two junction points x a and x b (Figure [2]), where a topo- 
logical phase transition occurs, we expect the emergence 
of finite energy ingap bound states localized in the vicin- 
ity of these points. However, if the distance d is taken 
to infinity, then the finite energy ingap states evolve into 
zero-energy Majorana bound states. 

First we shall examine the case where there is no 
phase mismatch along the junction and then turn to the 
general case. If we set all the phases equal to zero, 
we find that the energy eigenvalues satisfy the equa- 
tion E = || A| - \B\\e- K - (E)d . Notice that energy ap- 
pears in both sides of the equation. For a large dis- 
tance d, the energy tends to zero E — > and concomi- 
tantly K-(E) -> k_(0) = ||A| - \B\\/vH. By substi- 
tuting this approximate value for K- back to the equa- 
tion defining the energy, we directly obtain the expres- 
sion E = (|A| - |g|) e -(|AHB|)d/oft since | A | > 
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FIG. 2: (Color online) a. Topological-normal-topological 
junction, realized for the following choice of magnetic field 
and superconducting order parameter \Bi\ = \B\, \B m \ = |A], 
\B r \ = \B\, I Ai| = |A|, |A m | = \B\ and |A r | = |A|. b. 
Normal-topological-normal junction, realized for \Bi\ = |A|, 
\B m \ = \B\, \B r \ = |A|, |A,| = |S|, |A m | = |A| and 
|A, | = |S|. In both cases localized ingap states appear at 
the interface points x a and Xb where a topological quantum 
phase transition occurs. When the distance d becomes infi- 
nite, these localized states evolve into Majorana bound states 
7a and 76. 



We observe that the energy scales exponentially with the 
distance and it becomes zero for d — > 00, giving rise to 
the zero-energy Majorana bound states, one per interface 
point. In fact, for finite d, the finite energy splitting is 
a consequence of the hybridized Majorana bound states, 
that now have a finite overlap. 

In order to study the case where we keep the phases 
of the field and the order parameter intact, we shall as- 
sume that K±(E) ~ k±(0) = (|A| ± \B\)/vh. To obtain 
an approximate analytical expression for the energy of 
the bound states, we shall consider that the terms E, 
e -K+(o)d aiK j e -K-(o)d are £ same magnitude. In 

this manner, we may perform a perturbativc expansion 
up to second order in these three terms that will provide 
us with the following compact relation for the energy 



fl - <P, 



Jm cos 



J z cos I tpn 



(16) 



where we have defined the Josephson couplings Jm,z 
(| A| - \B\) (e~ K -(°) d ± e" K +W d )/2 and introduced 




B\)/[\A\ + \B 



cos x + cos UJ 



.(17) 



For a consistency check, we see that when the phases 
go to zero, T — > 1 and we obtain the anticipated result 
found earlier which yields E = (| A|-|£|)e~ (|AHB|)dM 



The term Jm cos 



-2 — . rJz cos (£i±£r — ^> m ) that ap- 
pears in our result, is in agreement with the findings 
of Ref.[53j]. The first coupling describes the usual 4n- 
periodic Josephson term and the second was recently 
highlighted by the aforementioned authors. We observe 
that our method not only retrieves the already estab- 
lished results but also additional information with sig- 
nificant physical consequences, which are encoded in the 
term T{ipi - <p m ,0l - fim)F{<Pr - <Pm>^r - &m)- Since 
this term includes the phases of the magnetic field "di !m ,r-, 
it provides the possibility of manipulating the Josephson 
response via the control of the magnetic field orientation 
along the junction. Based on reciprocity, we also pre- 
dict the generation of a spin current polarized along the 
z— direction, which can be controlled by tuning the su- 
perconducting phase difference along the junction. Both 
phenomena originate from the presence of the spin-orbit 
interaction of the semiconducting wire. 

To gain some more insight concerning these two phe- 
nomena, we shall consider the special case where ip m = 
Qpi + tp r )/2, i?m = (#i + #r)/2 and |A| >> \B\. In this 
case, everything depends on the phase differences ipi — ip r 
and §1 — *d r and we obtain 



, l|»l 
l ~2\K\ 



cos 



W - <Pr 



COS 



J' mC0S (^L\+J' 7 , 



(18) 



where J' M z correspond to the values that these 
quantities acquire when we take the limit |A| >> 
\B\. We readily observe that the Josephson cou- 
pling consists of three types of terms: cos 

and 



cos (<pi — <p r ) ^sincecos 2 [ Vl 2 yr ) 



l+cos(yj-y r ) 
~~ 2 



The first term is responsible 



cos cos 
for the 47r-periodic MBS Josephson effect, the second 
term describes a usual 27r-periodic Josephson effect and 
the last coupling describes a "magnetically" driven A-k- 
periodic Josephson effect or a superconducting phase 
driven 47r-periodic spin current. 



IV. NORMAL - TOPOLOGICAL - NORMAL 
(NTN) JUNCTION 

In this section wc shall consider a junction consisting 
of a normal-topological-normal sequence of phases. To 
model this type of junction we shall make the following 
choice for the order parameter and the magnetic field 
along the junction \B t \ = |A|, \B m \ = \B\, \B r \ = |A|, 
|A,| = \B\, \A m \ = |A| and |A r | = \B\. Once again, 
|A| and \B\ correspond to the values of the supercon- 
ducting gap and Zeeman field in the topological region. 
The methodology we use is identical to the one followed 
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in the previous section and we obtain the expression for 
the midgap state energy 



E = F{ipi - tpmj&i - ^ m )J r (v? r - ip m , d r - l?m) 



Jm cos 



Jz cos I - # r 



(19) 



where we have used the same definitions as previously. 
We observe that the results of this section can be related 
to the prior ones by just exchanging the phases <p -<-> '&. 
Again, if we take the special case <p m = (ipi + ip r )/2, 
#m = C#i +^, )/2 and the limit |A| >> we similarly 
obtain 



E = 



2IAI 



JmCos 



cos 



<Pl - ¥>r 



J' 7 . 



COS 



*I — *T 



In this case we 
kinds of couplings cos 



also retrieve 



(20) 

three different 
cos ($1 — i9 r ) and 



cos ( w 2^ r ) cos ( '2 r ). The first two terms are con- 
nected to 47r-periodic and 2-7r-periodic spin transport, 
while the novel type of "magnetic" Josephson coupling 
(third term) still appears in this case. 



V. MAGNETICALLY TUNED JOSEPHSON 
COUPLING AND ALTERNATIVE WAY OF 
BRAIDING 

Apart from the importance of exploring routes for en- 
gineering MBS and topological qubits, an equally signif- 
icant task is to invent new techiques for manipulating 
them and performing quantum computations. MBS offer 
due to their topological stability significant advantages 
as compared to the more conventional qubits based on 
either Josephson junctions [57j or spins in quantum dots 
[58| . However, for their manipulation new techniques 
have to be developed. One of the standard methods 
in the latter situation, is to create junctions (e.g. Y- 
junctions [H, EH) where a number of MBS are brought 
together, interact and afterwards are separated again in 
order for a read-out protocol to be implemented and yield 
the quantum computation output. Up to now, the only 
interaction considered to take place between two MBS 
was through the Josephson coupling 



Hm oc iji 7 r cos 



tpt - tp r 



(21) 



As a matter of fact, the superconducting phase of the 
order parameter has to be manipulated in order to braid 
the MBS. However, in this paper we revealed another 
type of Josephson coupling that exists in these systems 
and this is of the form 



H M " OC ljlj r COS 



tfl - lp r 



2 J 7 V 2 J ^ 
This novel Josephson coupling provides an additional de- 
gree of freedom to be manipulated, which is the orien- 
tation of the Zeeman magnetic field. In this manner, 
varying the superconducting order parameter phases is 
not necessary any more, since braiding can be performed 
by properly adjusting the magnetic field that each MBS 
feels. 



VI. CONCLUSIONS 

We presented a detailed analysis of the Josephson ef- 
fect in double junction sandwiches of topologically trivial 
and non-trivial phases of semiconductor-superconductor 
ID heterostructures. We demonstrated the existence of 
a novel type of Josepshon coupling that allows the ma- 
nipulation of the supercurrent via spatial variations of 
the Zeeman field along the junction. Conversely, a spin 
current generation is possible by varying the supercon- 
ducting order parameter phase differences. This addi- 
tional term opens a pathway for unambiguous Majorana 
bound state detection and at the same time constitutes 
an alternative platform for topological qubit operations. 

Note added: After we had already submitted this 
manuscript for the ICM2012 proceedings, two other pa- 
pers [59l |60| appeared, arriving independently to similar 
conclusions. 
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